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The low-energy monochromatic neutron emission from an 88 Y/Be source can be exploited to 
mimic the few keVnr nuclear recoils expected from low-mass Weakly Interacting Massive Particles 
(WIMPs) and coherent scattering of neutrinos off nuclei. Using this source, a <10% quenching factor 
is measured for sodium recoils below 24 keVnr in NaI[Tl]. This is considerably smaller than the 30% 
typically adopted in the interpretation of the DAM A/LIBRA dark matter experiment, resulting in 
a marked increase of its tension with other searches, under the standard set of phenomenological 
assumptions. The method is illustrated for other target materials (superheated and noble liquids). 

PACS numbers: 78.70.Ps, 29.40.Mc, 95.35.+d, 25.30.Pt 



Several recently observed anomalies in direct searches 
for dark matter PHI can be interpreted within the con- 
text of low-mass (few GeV/c 2 ) Weakly Interacting Mas- 
sive Particles (WIMPs). These hypothetical particles 
are expected to induce nuclear recoils [4] carrying just 
a few keV (keVnr) of energy in most targets. A recoil- 
ing nucleus produces a reduced amount of ionization or 
scintillation in detector materials, when compared to an 
electron recoil of the same energy. Careful study of this 
"quenching factor" is a necessary first step in the inter- 
pretation of dark matter detector data. 

A scarce number of calibration methods exist for this 
few keVnr nuclear recoil energy regime, and those typ- 
ically involve the use of accelerator or reactor facilities 
able to provide monochromatic neutron beams, and of 
small dedicated detectors. This Letter introduces a con- 
venient new method able to generate such recoils. The 
source is eminently fieldable, i.e., can be used in situ 
at large underground dark matter detectors, in this way 
incorporating specific instrumental effects into the mea- 
sured response. Planned devices aiming at detection of 
coherent neutrino scattering off target nuclei 0, Q can 
also benefit from its characteristics. 

Photo-neutron sources involve the gamma-induced dis- 
integration of a light nucleus with modest neutron bind- 
ing energy, such as 9 Be or 2 D [7|,l8j. Only a few relatively 
short-lived radioisotopes generate gammas above thresh- 
old for these reactions (1.67 MeV and 2.23 MeV, respec- 
tively) . The main advantage of these sources is the mono- 
energetic nature of the neutrons generated, typically in 
the few hundreds of keV. However, their neutron yield 
is several orders of magnitude lower than their gamma 
counterpart, due to the relatively small cross-section for 
the photo-disintegration reaction. On first impression, 
this large unbalance would complicate their use for cal- 
ibration of WIMP detectors, specially when considering 
that the ability to discriminate between gamma-induced 



electron recoils and neutron-induced nuclear recoils is of- 
ten very limited at the few keVnr expected. 
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FIG. 1. MCNP-PoliMi 10] simulated deposited energy by 
electron recoils (keVee) and nuclear recoils (keVnr) in a small 
NaI[Tl] scintillator irradiated by 88 Y/Be gammas and neu- 
trons, respectively. The 0.1 mCi source is 21 cm away from 
the detector, yielding an isotropic 400 neutrons/s and 7.4 x 10 6 
gammas/s. Several intermediate thicknesses of lead between 
source and detector are considered (see text). 

A way to bypass this issue is the addition of 15-20 cm 
of lead shielding around the source. Conveniently, this 
is a typical thickness already present around many large 
underground WIMP detectors. Lead combines optimal 
gamma shielding properties with one of the smallest neu- 
tron moderation powers in any material, due to its large 
atomic mass. Formally defined, this power is the prod- 
uct of neutron lethargy times the macroscopic scattering 
cross section, 2 x 10 -3 for lead and unity for the standard 
moderator, hydrogen Q. More intuitively, a ~100 keV 
neutron straggles in lead with a mean free path of ~3 cm 
between scatters, while losing just a maximum of 2% of 
its energy in each interaction, typically less. 

Fig. 1 exemplifies this approach for the case of a small 
NaI[Tl] scintillator in the presence of an 88 Y/Be source. 
The combination of 88 Y (T]/2= 107 d) and a beryllium 
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converter produces a dominant E n =152 keV neutron 
emission with just a ~0.5% E n =963 keV higher-energy 
neutron branch, from gamma rays at E 7 = 1.836 MeV and 
2.734 MeV, respectively The maximum recoil energy 
imparted to a sodium nucleus by the dominant neutrons 
is kinematically limited to (4 M m E n )/(M+m) 2 = 24 
keVnr, with M and m representing the sodium nucleus 
and neutron rest mass energies, respectively. The lead 
layer reduces the gamma flux arriving to the detector 
under test to manageable levels, while neutrons lose only 
a small fraction of their energy in traversing it. 
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sion from a 57 Co source (Fig. 2 inset), resulting in a light 
yield, integrated over 3 /is, in good agreement with that 



found in 11] via Compton scattering. These gammas 



are sufficiently penetrating to clear the 0.8 mm scintil- 
lator aluminum casing, while also avoiding skin effects 
from hydrated surface layers in the crystal |12| . 

Four pairs of runs were taken, with the 88 Y source sur- 
rounded alternatively by BeO or by a pillbox replica in 
aluminum metal, for a total of 90 hours in each configu- 
ration. Aluminum has gamma cross sections essentially 
identical to BeO for penetrating photons from this ra- 
dioisotope, but is inert for purposes of neutron emission. 
Two cuts were applied to the data during their analy- 
sis. First, an energy-dependent dead time was imposed 
following large energy depositions in the crystal, so as 
to avoid excess low-energy signals from the ensuing af- 
terglow (phosphorescence). Second, only events with a 
scintillation decay time r compatible with that from nu- 
clear and electron recoils were accepted. This range is 
50 ns < r < 400 ns for NaI[Tl] [12. Fig. 2 displays the 
cumulative residual energy spectrum obtained by taking 
the difference between 88 Y/Be (neutron plus gamma) and 
88 Y/A1 (gamma only) runs. The low-energy excess from 
neutron-induced nuclear recoils is readily visible in this 
residual, and present for all pairs of runs. 



FIG. 2. Measured response of NaI[Tl] to low-energy nu- 
clear recoils from an 88 Y/Be source (see text). A solid his- 
togram is the predicted response, allowing no free parame- 
ters, obtained by adopting the quenching factors for sodium 
and iodine recoils recently measured using 2.2 MeV neutron 
scattering from a D-D generator (a monotonically decreas- 
ing Qnu with decreasing recoil energy, and Qi — 0.04 [Ilj]). 
A dashed histogram, in large disagreement with present data, 
employs quenching factors typically used in the interpretation 
of DAMA/LIBRA results (Q Na = 0.3, Q/ = 0.09 [ljj). 

Fig. 2 demonstrates the experimental application of 
this concept. A small (1.9 cm diameter, 5.1 cm 
long) NaI[Tl] scintillator from Proteus/ Amcrys was sur- 
rounded by 5 cm of lead, except for on one side of a 
rectangular box enclosure, 20 cm thick. Outside of this 
thicker lead wall, a small cylindrical "pillbox" (4.8 cm 
diameter, 2.5 cm tall) of BeO obtained from American 
Beryllia was used to encapsulate an 88 Y button source, 
procured from Eckert & Ziegler Isotope Products. The 
activity of this source at the time of the measurement was 
approximately 50 /iCi, with an isotropic neutron yield of 
~330 n/s, characterized with a 3 He counter surrounded 
by polyethylene moderator, the response of which was 
calculated via MCNP simulation. 

The same data acquisition (DAQ) system as in 
was employed, set to trigger on single photoelectrons 
(PE) from a Hamamatsu ultra-bialkali photomultiplier 
directly coupled to the scintillator. The rapid triggering 
rate (~800 Hz in the presence of the source) was well- 
below the maximum throughput of this DAQ. An energy 
scale was established using the 14.4 keV gamma emis- 
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FIG. 3. DAMA/LIBRA ROI in spin-independent coupling 
vs. WIMP mass phase space, resulting from the adoption of 
the quenching factors derived from present measurements and 
those in [TTJ] , compared to the ROI generated by a conven- 
tional Q Na = 0.3, Q/ = 0.09 0. Solid contours correspond to 
a 90% C.L. region, dashed to 99% C.L. A Maxwellian galac- 
tic halo is assumed, with local parameters vo =230 km/s, 
v esc =600 km/s, p =0.3 GeV/c 2 cm 3 . A CoGeNT ROI is 
from 14] and includes a correction for residual surface back- 
grounds as in [15]. A possible ROI for CDMS is from p]|. 
ROIs and limits from CRESST 0, CDMS [U, EDELWEISS 
[3, TEXONO [3, and XENON100 [H are also shown. 
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A solid histogram in Fig. 2 represents the simulated 
nuclear recoil response following adoption of the NaI[Tl] 
quenching factors most recently measured in [ll[ , also in- 
troducing a simple model for the energy resolution based 
on the measured response to 57 Co. The simulation in- 
cludes self-shielding at the source and the (negligible) ef- 
fect of known chemical impurities in the lead used. This 
prediction is found to be in good agreement with the 
data, even before any attempt at optimizing the energy 
dependence of the quenching factor or the resolution. A 
dashed histogram represents the prediction using the con- 
siderably larger quenching factors [13] typically adopted 
in a WIMP interpretation of the DAM A/LIBRA results. 
It is at clear variance with present data. It should be 
noted that these conventional values, QNa = 0.30 ±0.01 
and Qi = 0.09 ± 0.01 arise from an averaged response 
over 6.5 to 97 keVnr and 22 to 330 keVnr, respectively 
[l3| , i.e., from recoil energy regions much broader than 
those of actual interest for low-mass WIMP studies. 
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FIG. 4. Measured response of superheated CF 3 I to 241 Am/Be 
(<E n >= 4.5 MeV) and 88 Y/Be neutrons (E n = 0.152 MeV). 
As expected upon removal of the Be converter, only a minimal 
sensitivity to 88 Y gammas is observed, starting at very high 
degrees of superheat [22|, HH- Inset: simulated distribution of 
88 Y/Be nuclear recoil energies explored by the decompression 
protocol. These distributions can be skewed towards higher 
or lower mean energies by adjusting the decompression rate 
of the bubble chamber. No contribution from iodine recoils is 
expected in the slow decompression protocol employed. 



Fig. 3 shows the sizable displacement in the 
DAM A/LIBRA region of interest (ROI) that stems from 
the adoption of smaller quenching factors like those mea- 
sured here and in [11]. The best-fit WIMP mass able to 
reproduce the DAM A/LIBRA annual modulation effect 
is shifted from ~ 10 GeV/c 2 to ~ 18 GeV/c 2 (crosses in 
the figure), pushing this ROI further into the region of pa- 
rameter space excluded by numerous other experiments, 
and away from other anomalies. The reader is referred 
to [ll| for an additional discussion on astrophysical un- 
certainties perhaps capable of easing this now enhanced 



tension. Systematic effects affecting previous measure- 
ments of quenching factors at few keVnr, able to generate 
artificially large values, are described in (ill l2lj|. 

Fig. 4 displays results from the 88 Y/Be irradiation of 
a small bubble chamber built by the author to study the 
response of moderately superheated CF3I to low-energy 
nuclear recoils. In this particular case, no lead shielding 
is necessary, given the very large intrinsic insensitivity of 
this target to electron recoils [22]. The protocol employed 
in these measurements is that of a slow controlled decom- 
pression of the chamber in the presence of the source, 
with or without Be converter, effectively scanning the 
theoretically-predicted nuclear recoil energy threshold. 
As this threshold is lowered, a "hot border" [23| is en- 
countered when the maximum recoil energy produced by 
the monochromatic neutrons is reached, inducing bubble 
nucleations. The pressure cycle is then repeated. This 
protocol allows to extract the exact energy-dependence 
for this threshold from the analysis of the data [24j . 

Lastly, Fig. 5 shows the expected response to a lead- 
shielded 88 Y/Be source in a small liquid xenon (LXe) 
calibration cell comparable to that in [25| . Overlapped is 
the very similar expected response to low-mass WIMPs 
like those of interest vis-a-vis a dark matter interpreta- 
tion of the anomalies in The use °f this source 
situ with a larger LXe chamber such as the XENON100 
dark matter detector would incorporate any instrumental 
effects (e.g., light collection efficiency, hardware trigger 
configuration, etc.) into the calibration of the detector 
response, but would also lead to higher-energy signals 
from multiple neutron scattering. These can however 
be identified through the good ~1 mm spatial recon- 
struction of this detector [26]. This straightforward cal- 
ibration would provide a definitive test of the response 
of the XENON100 apparatus to low-mass WIMPs, re- 
solving the controversy around its claimed sensitivity to 
WIMP masses below - 12 GeV/c 2 
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FIG. 5. Simulated response of a small LXe detector to an 
88 Y/Be source (nuclear recoil energies in keVnr, electron re- 
coil energies from gamma component in keVee). A very simi- 
lar spectral response to low-mass WIMPs like those of present 
interest is shown, normalized to an arbitrary cross-section. 
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In conclusion, a new and convenient calibration 
method able to replicate the expected signals from low- 
mass WIMP scattering has been demonstrated. Its use 
with a NaI[Tl] target confirms recent quenching factor 
measurements by the author using an independent ap- 
proach [ll[ . These results substantially increase the diffi- 
culty in interpreting the DAM A/LIBRA anomaly within 
a WIMP dark matter context, under the standard par- 
ticle physics and astrophysics phenomenological conven- 



tions. In an upcoming publication 88 Y/Be will be 
used to study the decay time of scintillation light from 
nuclear recoils in CsI[Na], for the few keVnr energy range 
expected from coherent neutrino scattering off nuclei, an- 
other possible realm of application. 
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